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After the start of the LHC, a plethora of novel observables for jet tomography in heavy-ion colli- 
sions has appeared. Many of these studies have yielded counter-intuitive null results of apparently 
unmodified jets, which have sparked (sometimes exotic) theoretical efforts to explain these findings. 
However, it has to be realized that almost all current high Pt observables measure conditional proba- 
bilities of events, not probabilities. Thus, the correct starting point for their theoretical understand- 
ing is Bayes' formula, and the biases introduced by the conditioning are crucial to understanding the 
outcome. Once this is introduced properly into the modelling process, the counter-intuitive results 
are seen to find a natural explanation in terms of various biases and the puzzles largely disappear. 
In this work, a conceptual framework to classify the various observables according to the types of 
bias they are sensitive to is presented and illustrated with a large number of case studies ranging 
from simple jet finding to 2-1-1 dihadron triggered correlations. 



PACS numbers: 25.75.-q,25.75.Gz 



I. INTRODUCTION 



The idea of doing jet tomography in ultrarelativistic 
heavy-ion (A-A) colhsions, i.e. to utihze hard processes 
taking place along the creation of a soft bulk medium 
to probe both the medium geometry and its degrees of 
freedom has been proposed many years ago [ll-@| . At the 
Brookhaven Relativistic Heavy Ion Collider (RHIC), the 
observables considered to probe this physics was initially 
the nuclear suppression factor of single inclusive hadrons 
Raa 0, 111 and the suppression factor Iaa of hard back- 
to-back dihadron correlations [1, [l0| • 

Recent high statistics runs at RHIC as well as the 
significantly larger kinematic reach of heavy-ion exper- 
iments at the CERN Large Hadron Collider (LHC) have 
led to a sometimes bewildering plethora of new high Pt 
observables, in particular also observables involving jet 
reconstruction, among them dijet imbalance measure- 
ments [ill [T2|. iet-hadron (jet-h) correlation s fisl . h-jet 
correlations , jet fragmentation function [1^ and jet 
shapes [l^ or observales utilizing rare electroweak trig- 
gers such as 7-h correlations [13: [3 or 7-jet correlations 

To add to the complexity, jet definitions vary from 
calorimetric jets which are not corrected for background 
fluctuations as used e.g. in [ll| to combined track/tower 
jets with Pt cuts imposed on constituents and a hard 
tower trigger condition imposed as used e.g. in 13]. 

In this situation, it is fairly difficult to assemble a pic- 
ture of what information the various observables actu- 
ally carry, to what degree they are mutually consistent 
and what features of models they constrain. The aim 
of this paper is to improve on this situation by provid- 
ing a clear conceptual framework in which similarities 
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and differences between the various observables become 
transparent. 

The key observation for this is that the vast majority 
of observables (with the exception of nuclear modifica- 
tion factors) are measurements of a conditional proba- 
bility given a trigger condition. The fundamental reason 
for this is that both hard and electroweak processes are 
rare, i.e. if there would be no selection of the subclass of 
events containing hard processes, the background of soft 
bulk medium physics would dilute all signatures of hard 
probes to the point that they would no longer be observ- 
able. However, conditional probabilities are well known 
to be frequently non-intuitive, and the natural starting 
point for analyzing them is Bayes' formula, which will be 
utilized in the following. 



II. OBSERVABLES AND CONDITIONAL 
PROBABILITIES 

A. General considerations 

In perturbative Quantum Chromodynamics (pQCD), 
the rate of hard scattering processes can be computed 
with reasonable accuracy once the momentum transfer 
in the scattering process exceeds a few GeV. The uncer- 
tainty principle allows to estimate the timescale for the 
hard reaction as r ~ iS/Q^ where E is the energy scale of 
the final state partons and Q ^ 0{E) the virtuality scale. 
Inserting typical numbers, one finds that hard processes 
occur before a soft medium can be formed, which is the 
reason that the pQCD computation of hard processes can 
safely be assumed to factorize from any medium physics. 
This property makes high Pt observables a meaningful 
tomographic probe. 

The highly virtual back-to-back partons subsequently 
undergo a final state shower evolution in which the virtu- 
ality scale decreases from its initial high value to a non- 
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perturbative scale via the branching into additional par- 
tons. This process in vacuum is well described by MC for- 
mulations such as the PYSHOW algorithm of PYTHIA 
[20j . Once at the non-perturbative scale, the parton 
shower hadronizes and becomes a coUimated spray of 
hadrons. Jet clustering algorithms such as anti-fcr or 
SIS-cone as provided e.g. by the FastJet package [2l| 
aim to 'undo' the QCD shower evolution and turn the 
spray of hadrons again into a 'jet', i.e. an object which is 
a reasonable proxy for the original parton largely free of 
the complications of shower evolution and hadronization 
and sensitive to hard physics only. 

Measurements of hard probes in the context of heavy 
ion collisions aim at answering the question how the 
medium modifies this evolution, i.e. in what way the 
properties of the shower are different if it evolves inside 
a medium. If a jet contains n hadrons, since the posi- 
tion space information can not be resolved, the complete 
experimentally available information about the jet is con- 
tained in the momentum space density pn[Pi, P2, ■■■Pn) 
and in the knowledge of hadron identities. However, cur- 
rently the focus is on measurements of the single parti- 
cle distribution pi(Fi) = / dP2 . . . dPnPn{Pi,P2, ■ ■ - Pn), 
usually represented as longitudinal and transverse mo- 
mentum spectra of particles with respect to the jet axis. 
In the future measurements may also include intra-jet 
correlations as given e.g. by two particle correlation 
C2{Pi,P2) and three particle correlations C3(Pi,P2,^3) 
or in terms of subjet fractions. This information may be 
represented in different form, for instance the integrated 
jet shape 
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(the integrated flux of energy as a function of angle r 
with the jet axis of a jet of radius R, normalized to the 
total jet energy) is computable from the angular distribu- 
tion of hadrons at given energy dN/ d(j)dE (as for instance 
obtained from a correlation measurement) as 
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(2) 



It is thus sufficient to measure one representation of 
the single particle distribution, different representations 
contain redundant information. 

If we use the notation that P{A\B) stands for the 
probability of event A occurring given another event B, 
the computation of the probability of observing shower 
properties S (for instance the probability of measuring 
a shower hadron between 2 and 3 GeV) given a set of 
trigger conditions T (for example given that a jet is clus- 
tered in an energy between 100 and 150 GeV) is written 
as P{S\T, M) where M stands for the particular model 
in which the calculation is carried out. B ayes' formula 
then allows to compute this as 
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In words, the probability for observing shower proper- 
ties S given a trigger T is the product of the probability 
to fulfill the trigger condition in a shower with property 
S times the probability to generate a shower 5, divided 
by the probability to generate a trigger independent if 
property S is realized or not. Since a rate is obtained 
by multiplying a probability with a repetition frequency, 
the whole language trivially generalizes to event rates or 
particle spectra. 

What is measured is usually the left hand side of the 
equation, somtimes also the denominator of the right 
hand side (which corresponds to the rate at which the 
trigger condition is fulfilled). Eq. ^ states then that in 
a large class of measurements, the medium modification 
as computable in a model P{S\M) is not be observed 
directly, but rather is distorted through a bias factor 
^]^('r\ 'M) '^hich is characterized by the trigger condition 
T . This bias can vary a lot, for instance the requirement 
to find a 100 GeV calorimetric jet leads to a very differ- 
ent bias than the requirement to find a 20 GeV charged 
hadron, but as these examples indicate, the formalism ap- 
plies as well to jet finding followed by an analysis of the 
fragmentation pattern of the clustered jet psi] in which 
case the jet finding constitutes the trigger condition and 
the observable is the longitudinal momentum spectrum 
of the shower, as to Iaa in triggered h-h correlations in 
which the requirement to find a hard hadron constitutes 
the trigger condition and the ratio of longitudinal mo- 
mentum spectra of correlated hadrons in medium over 
vacuum is the observable. 

This suggests a clear strategy to make the informa- 
tion content of measurements apparent and comparable: 
Measure the observable (e.g. the single particle distribu- 
tion of jet constituents) in the same representation in all 
measurements and view the different trigger condition as 
a variation of the bias factor. Tomographic information 
is then contained in the way the observable responds to 
a change of the bias factor. 



B. Theoretical formulation of in- medium showers 

As discussed in detail in '22'] , modelling of the medium 
modification of a shower involves a procedure to compute 
the medium-modified fragmentation function (MMFF). 
The MMFF can be written in the rather general form 
A^/.(z,£^,gg|Ti(C),r2(C),...r„(C)), where it describes 
the distribution of hadrons h given a parton i with ini- 
tial energy E and initial virtuality Qg where the hadron 
energy Eh = zE and the parton has traversed a medium 
along the path C, where Ti{C,) are the medium transport 
coefficients relevant for the process. 

Since the MMFF should approach the usual vacuum 
fragmentation function when the transport coefficients 
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vanish, the properties of a vacuum shower are largely 
determined by just three parameters — the shower- 
initiating parton type i, its initial energy E and vir- 
tuality Qq. In contrast, the determination of medium 
modifications in principles require n different functions 
Ti{(). However, it turns out that in practice three are 
most relevant: q (the virtuality growth per unit path- 
length), e (the energy loss per unit pathlength) and 
62 (the variance of the energy loss) [23|. Moreover, in 
many models it turns out that the full functional de- 
pendence of the transport coefficient is not needed but 
rather the first one or two moments Mi = J d(Ti{() 
and M2 = J d((Ti{() are to good accuracy sufficient 
[U, [2^ . This implies that the medium modification of 
a shower can be characterized reasonably well by the 
set of Ml (g) , M2 (q) , Ml (e) , M2 (e) , Mi (62) , M2 (ea ) which 
now contain all tomographic information. 

Thus, ideally one would like to compare 
Di^ii{z, E, QQ\Mi{q), . . . ) with a measurement to 
deduce the tomographic information on the properties 
of the medium. 

A number of models for the computation of the MMFF 
are proposed. Historically, the computation has often 
been based on the leading parton energy loss approxima- 
tion in which the virtuality evolution of the shower is not 
treated explicitly and the focus is only on induced radi- 
ation from the leading parton 0, 0, [26l - l29| . Since this 
approximation is not well suited for the interpretation in 
terms of conditional probabilities, we will not consider 
it here. Alternatively, Monte-Carlo (MC) codes for in- 
medium shower evolution [2^ Isol ^ 
as well as analytical approaches 



Consider the term P{T\M) which can be written as 
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3c I , parton cascade [34 
35, exist. 



C. Initial state and final state biases 

In order to directly test a model of jet quenching, it 
would be desirable if an observable could be constructed 
in such a way that the vacuum shower model parame- 
ters {i,E,Qo) or the medium parameter moments take 
fixed values. In this case, the theoretical model would 
only ever need to consider events which fulfill the trig- 
ger condition by construction, rendering the bias factor 
^(t'^a/)"' identically unity, which simplifies the compu- 
tation tremendously. This is the reason schematic inves- 
tigations and toy models follow this strategy. In other 
words, if one could prepare a situation in which a quark 
with specified energy propagates through a given length 
of medium with given density, jet tomography through 
comparison of experiment and theory would be easy to 
do. 

Unfortunately, experimental measurements are hardly 
ever conditioning on initial state properties of the shower, 
in which case the bias factor is different from unity and 
a model to compute the MMFF is insufficient to com- 
pare with data. Instead, experimental trigger conditions 
usually key on some property of the observed final state 
after shower evolution and hadronization. 



using the fact that probabilities normalize to unity. 
Eq. (lU states that in order to compute the rate at which 
the trigger condition is fulfilled, we need not only com- 
pute the shower S exhibiting a particular property we're 
interested in but in fact all possible shower configurations 
and medium modifications S' which are allowed by the 
physics of the collision and do an appropriate sum over 
them. It is this need to compute all possible initial con- 
figurations and check them for the trigger condition in 
the final state which makes a proper computation vastly 
more complicated than a toy model estimate. 

In practical terms, this means that in order to compute 
observables which can be compared with experiment, an 
in-medium shower model needs to be embedded into a 
framework simulating the hard process and the evolution 
of^he surrounding medium (for a detailed discussion see 

The final state trigger condition than maps (in a 
model- and embedding-dependent way) into distributions 
in the space of initial shower parameters. 



D. Monte-Carlo treatment of jet quenching 

Let us for illustration consider a MC description of jet 
quenching. Biases are taken into account by generating 
events according to the full abailable space of initial pa- 
rameters with correct weight assigned to the individual 
contributions, then searching which of these events fulfill 
the trigger condition in the final state and analyzing only 
this subset of event to obtain the observable. Pictorially 
this is shown in Fig. [TJ Ultimately interesting for the ob- 
servable are only the two shaded regions, i.e. the class of 
events which fulfills the trigger T and the class of events 
which fulfills T and shows property S. 




all possible initial parameters 



FIG. 1: A schematic illustration of initial parameter space 
sampling for a conditional probability observable. 

However, the computational problem is that the full 
range of events generated by sampling all the available 
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initial parameter space is usually so huge that a naive ap- 
plication of the above strategy is bound to be so slow that 
it is useless in practice. The challenge resulting from this 
is to introduce an intemediate layer, i.e. to understand 
the bias structure in such a way that only initial param- 
eter ranges are sampled which have a reasonable chance 
to lead to a trigger in the final state. Pictorially, this 
corresponds to drawing the dashed line as closely as pos- 
sible to the intersecting circles without actually cutting 
parameter space out of a circle (which would introduce 
an unphysical sampling bias). In this way, computations 
become feasible. This illustrates that good knowledge 
of the mapping of final state conditions to initial state 
parameters in terms of biases is not only conceptually 
important, but also has consequences of immediate prac- 
tical value. 



III. TYPES OF BIASES 

Following the discussion in [s^], we can classify the 
various biases induced by a trigger condition on the final 
state of a hard event as follows: First, there are biases 
on the structure of the hard pQCD event itself which act 
even in vacuum. These have to do with the relation be- 
tween hadronic (or jet) and parton kinematics dependent 
on parton type. Once a medium is present, the correla- 
tion of the strength of the medium modification with the 
density of the medium and the time spent in the medium 
leads to additional biases on the reaction geometry. Since 
all these biases act on the hard event itself rather than 
the final state shower, they affect both trigger side and 
away side simultaneously. This can be contrasted with 
shower biases, which affect the structure of the shower 
evolution itself and do not bias the kinematics or posi- 
tion of the hard event and are thus always only relevant 
for the trigger side. 

In this section, we review qualitatively the effects of 
the most relevant biases, which we illustrate later with 
case studies in a full modelling framework. 



A. Biases in vacuum showers 

Neither a hadron nor a jet typically contain all the 
initial parton energy E. In the case of a hadron, this 
is because of the production of subleading hadrons as 
well as hadron species which are not registered by the 
detector in the shower. In the case of a jet, the reason is 
typically the production of hadrons at large angles with 
the jet axis which correspond to energy flow outside the 
jet radius R, but for instance in charged jets also neutral 
hadron production in the shower constitute an energy 
component not part of the jet. 

For both jet and hadron, the relation of observed 
energy to parton energy can be written into the form 
Eobs — ^had/jetE/- Typically, the chief difference between 
jet and hadron observation is that a jet tends to recover 



a higher fraction of the parton energy than a single hard 
hadron, i.e. (zjet) > (zhad) where the average is done 
over many showers with a fixed parton energy E. 

This is illustrated in Fig. [5] where P{z), the prob- 
ability to observe the fraction z of the original en- 
ergy of a 20 GeV quark in the final state is shown for 
three different objects: 1) the leading hadron if it is 
7r+,7r-,7r",K+,K-,p or p 2) a STAR jet definition [l3] 
where all particles which are ,Tr~ ,Tr'^ , , ,p or p 
or 7 and have Pt > 2 GeV are clustered usign the anti- 
kx algorithm with a radius of R = 0.4 and 3) an ideal 
jet definition where all particles, regardless of FID or Pt, 
are clustered with anti-Zcy using R — 0.4. 

20 GeV quark 




z 



FIG. 2: The probability density P{z) to observe a trigger 
object with fraction z = Eobs/E given an initial parton energy 
E and an observed trigger energy Eobs for various possible 
trigger objects, shown for the example of a fragmenting 20 
GeV quark. 

It is evident that the leading hadron in this kinematical 
regime typically carries only about 15% of the original 
parton energy whereas on the other end of the spectrum 
clustering into a jet ideally recovers typically 95% of the 
energy. Jet definitions matching realistic experimental 
conditions fall between the two cases. 

A kinematic bias arises then because in an experimen- 
tal context P{z) is typically not probed for fixed parton 
energy, but rather folded with the steeply falling primary 
parton production spectrum which can be computed in 
pQCD and typically falls approximately like a power 
1/Pt with n — 7.. 8 at RHIC kinematics and n ~ 4.. 5 
at the LHC. A trigger energy requirement then demands 
a fixed Eobs = zE where both z and E are allowed to 
vary event by event. For the ideal jet described above 
where P{z) 6{z — 1), the bias is negligible and Eobs 
approximately corresponds to the parton energy. For a 
hadron trigger however, both E and z prefer to be in- 
dividually small, yet their product is forced to a certain 
value. As a result, Eobs maps to a characteristic range in 
E which depends on n and the details of P{z), i.e. the 
distribution of parton energies contributing to a trigger 
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is no longer the primary pQCD spectrum but becomes 
biased. 

Another part of the kinematic bias is related to the fact 
that due to higher order pQCD effects and nuclear initial 
state effects a hard parton pair is never exactly back to 
back. These effects can be approximated by introducing 
a randomly oriented vector kt with a Gaussian distribu- 
tion in magnitude which is added to the pair momenta. A 
trigger condition then biases this a priori randomly ori- 
ented vector to be pointing towards the trigger direction 

The parton type bias then has to do with the fact 
that the functional form of P{z) depends on the shower- 
initiating parton type: On average, quarks fragment into 
harder and more coUimated showers than gluons. As a re- 
sult, any trigger condition corresponding to an observed 
energy is more likely to be fulfilled by a quark than by 
a gluon. Thus, on the trigger side the fraction of quark 
jets is generically enhanced as compared to an unbiased 
pQCD spectrum. How the bias acts on the away side 
depends on the kinematical situation. In a regime where 
the subprocess qg — > qg is dominant, enhancing the near 
side quark fraction biases the away side towards gluon 
jets [33 which is relevant for instance for the 5-20 GeV 
momentum regime at RHIC. 



B. Biases in tlie medium 

Medium modifications to the shower structure generi- 
cally tend to equilibrate the shower, i.e. they drive the 
kinematical properties of shower partons closer to those 
of medium partons. This implies that medium-modified 
showers are softer and broader, i.e. more weight in P{z) 
shifts to lower z. As a result, the kinematical bias is 
changed in a medium, the same Eots maps on average to 
a higher E in a medium than in a vacuum. 

The strength of the medium modification is (up to co- 
herence effects which are important in detail) driven by 
the number of interactions with the medium, which is a 
function of the medium density, the coupling strength of 
shower partons to the medium and the time/length of 
the shower spent in medium. Out of these, the coupling 
strength is relevant for a modification of the parton type 
bias by the medium: Since gluons interact with a factor of 
9/4 more strongly with color charges, the medium mod- 
ification of gluon jets is correspondingly stronger than 
that for quark jets. For this reason, triggered objects are 
even more biased to be quark jets than this is already 
the case in vacuum. This effect is sometimes referred to 
as gluon filtering. 

The combined effect of medium density and pathlength 
of a parton through the medium leads to a geometrical 
bias on the position of the vertex leading to the triggered 
event in the transverse plane in position space. Vertices 
leading to triggered events have a tendency to be close 
to the medium surface, with the trigger parton travel- 
ling outward. This implies that the same effect biases 



the away side parton to have a longer than average path- 
length in the medium. 



C. Shower biases 

While all biases discussed so far affect properties of the 
hard event itself, and thus refer equally to near and away 
side, there are also biases which affect the trigger parton 
side only. Those are here referred to as shower biases. 
For instance, requiring that a single hard hadron is pro- 
duced in a shower restricts the phase space for associate 
hadron production via the conservation of energy and 
momentum. Generically, shower biases make observables 
more robust against medium modifications, as a shower 
bias implies that there are properties of the shower which 
are by the trigger condition protected against medium 
modifications. 



IV. MODEL DESCRIPTION 

In order to illustrate the qualitative remarks made 
above quantitatively, we will in the following show results 
obtained with the in-medium shower evolution code Ya- 
JEM [H [13] in its latest version YaJEM-DE 38] which 
gives a fair account of a large number of observables both 
at RHIC and at LHC 
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YaJEM is a tool to obtain the MMFF given initial 
parton energy and a path through the medium, hence 
for a complete model description of a hard process in a 
medium also the medium evolution and the pQCD pro- 
cess have to be taken into account. 



A. The perturbative hard process 

Any simulation of hard events inside a heavy-ion colli- 
sion must start with the computation of the hard process 
itself. In LO pQCD, the production of two hard partons 
fc, / is described by 



da 



AB^kl+X 



dp^dyidy2 



^Xifi/A{xi,Q'^)x2fj/B(x2,Q'^ 



dt 



(5) 

where A and B stand for the colliding objects (protons 
or nuclei) and yi{2) is the rapidity of parton k{l). The 
distribution function of a parton type i in A at a mo- 
mentum fraction xi and a factorization scale Q ^ pT m 
.fi/Ai^i^Q'^)- The distribution functions are different for 
free protons [s^, and nucleons in nuclei [ill - lisj . The 
fractional momenta of the colliding partons i, j are given 
by a;i 2 = ^ (exp[±yi] -I- exp[±?;2])- Expressions for the 

pQCD subprocesses ^^^-^ — {s,i,u) as a function of the 
parton Mandelstam variables s, t and it can be found e.g. 
in H. 
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To account for various effects, including liigher order 
pQCD radiation, transverse motion of partons in the nu- 
cleon (nuclear) wave function and effectively also the fact 
that hadronization is not a collinear process, the distri- 
bution is commonly folded with an intrinsic transverse 
momentum kr with a Gaussian distribution, thus creat- 
ing a momentum imbalance between the two partons as 

PTi + PT2 = kx- 

In a MC description of the process, Eq. ([SJ is sampled 
to generate the parton type and momentum of the back- 
to-back pair. Subsequently the intrinsic kx imbalance 
is sampled and added to the parton pair momentum. In 
correlation studies, one of the partons is randomly picked 
as a trigger candidate. 



B. Medium-modified fragmentation 

Hard vertices are assumed to be distributed with a bi- 
nary overlap profile as appropriate for LO pQCD parton 
production, i.e. the a priori probability density for find- 
ing a vertex in the transverse (x, y) plane is given 



and 



r^(ro + b/2)r^(ro ~b/2) 



(6) 



TAA{h) 

where the thickness function is given in terms of Woods- 
Saxon distributions of the the nuclear density p^(r, z) 
as Ta{y) = J dzpA{r,z) and TAA(i>) is the standard nu- 
clear overlap function Taa (b) — J (Ps Ta{s)Ta (s — b) for 
impact parameter b. 

In the MC procedure, we place the parton pair at a 
probabilistically sampled vertex {xq^i/q) sampled from 
this distribution with a random orientation (jj with re- 
spect to the reaction plane. In this work rotate the event 
for the purpose of extracting vertex distributions such 
that the vector of the trigger candidate parton defines 
the —X direction, in studies of explict dependence of ob- 
servables on the angle of the parton with the various w„ 
event planes we would use the relevant event plane angle 
instead and would consider only parton propagation with 
a set angle to the event plane. 

The event is now embedded into a hydrodynamical de- 
scription of the medium ( [i^ for the RHIC case and the 
extrapolation of this scenario to larger y/s in the LHC 
case [52|) which allows to extract e.g. the energy density 
e(C) at any point of the propagating parton path 

In the absence of a medium, YaJEM is identical to 
the PYSHOW algorithm which evolves partons as 
a series of a — )• 6c branchings in the energy fraction 
z = Eb/Ea and the virtuality t = \n{Ql)/ Aq(jjj with 
Aqcd = 0(300) MeV. In YaJEM, it is assumed that the 
virtuality and energy Ea of any intermediate shower 
parton a is modified by the medium via two transport 
coeffients, q and e as 



dcm 



(7) 



AE„ 



dCe(C). 



(8) 



To evaluate these equations requires a mapping of the 
shower evolution of PYSHOW in momentum space to the 
hydrodynamical evolution in position space and a model 
of the transport coefficients as a function of thermody- 
namical properties of the medium. 

The temporal structure of the shower evolution can be 
parametrically recovered by uncertainty arguments. The 
mean lifetime of a virtual parton b coming from a parent 
a is hence given as 



Eb Eb 



(9) 



In the MC simulation of the shower, the actual life- 
time is determined from this mean value according to 
the probability distribution 



P{Tb) = exp 



Tb 

in) 



(10) 



For the relation between transport coefficients and hy- 
drodynamical parameters. 



m (C) = Kq [Ke] ■ 2 • [e(C)]3/^(cosh piO - sinh p(C) cos ^) 

. . ^^^^ 

is assumed where p is the transverse flow rapidity of 
the medium, -0 the angle between parton direction and 
medium flow direction and Kq and Ke are two free pa- 
rameters parametrizing the strength of the coupling of 
medium and shower partons. Following the procedure in 
fssj . these are adjusted to Kq — 0.8K and Ke — O.IK 
(corresponding to about a 10% elastic energy loss contri- 
bution leading to direct energy transfer into the medium) 
and K is fit to the nuclear suppression factor Raa in 0- 
10% central Au-Au collisions at RHIC. 

Note that in this work, any results presented are in- 
tended as illustration of qualitative effects and the order 
of magnitude of various biases. Hence no attempt is made 
to obtain a good fit to any other data set by either fitting 
K to a more extended set of data or by exploiting the 
freedom to choose a different fluid dynamical description 
of the medium. For this reason comparison with data 
(where it is available) is left for future work. 

Changing the kinematics of evolving shower partons 
according to Eqs. (|7l8p in YaJEM results in a medium- 
modified parton shower. The resulting distribution of 
partons is then passed to the Lund model (46j to compute 
the non-perturbative hadronization. 
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C. Analysis 

The resulting output in terms of medium-modified 
hadron showers is now analyzed if the trigger condition is 
fulfilled. Since the event record at this point contains the 
full information on hadron PID and momenta, in prin- 
ciple any set of cuts can be evaluated (in practice the 
statistics may become too low). 

In case the trigger is a hadron, the test for the trig- 
ger condition is trivial. In the case of a jet trigger, the 
resulting event record is clustered with the anti-fcr algo- 
rithm of the Fast Jet package [2l|. At this point parti- 
cles computed from a bulk hydrodynamical event could 
be inserted and clustered together with the hard event 
to study the influence of background fluctuations. This 
however is computationally very expensive and not done 
here — throughout this work, it is assumed that any 
background fluctuations are sufficiently trivial to be re- 
moved. Dependent on the trigger conditions, this niay 
not be a good assumption in practice (see e.g. [U Hsj 
for studies of the influence of the soft background on ob- 
servables) . 

In the case a jet trigger in combination with conditions 
on PID or constituent momenta is required, all particles 
not fulfilling the conditions are removed before clustering 
is done. 



D. Relevance of the results 

Despite the fact that the following case studies are per- 
formed with a specific parton-medium interaction model, 
YaJEM-DE, and for a specific choice of medium evolu- 
tion, the qualitative conclusions drawn about the role 
of biases in hard observables require a significantly less 
stringent set of assumptions. 

The medium-induced biases will appear acting into 
the same direction as illustrated by YaJEM-DE for any 
model which has the following characteristics: 1) the 
medium on average softens fragmentation in a shower 
2) the medium on average broadens the transverse dis- 
tribution of hadrons in a shower 3) the effects of soften- 
ing and broadening increase monotonously with medium 
density and in-medium pathlength 4) gluons couple more 
strongly to the medium than quarks. 

Most current models of parton-medium interactions 
exhibit these traits, so the following qualitative state- 
ments can be expected to hold fairly generically. How- 
ever, quantitatively the relative strength of different bi- 
ases depends on specific model assumptions. 

V. SHOWER BIASES 

A trigger condition may refer to one or both showers 
generated in a back-to-back hard event. If the shower on 
the trigger side is studied, in general the trigger condi- 
tion biases the shower evolution itself. The shower bias 



is however absent when the trigger condition is evalu- 
ated for one parton (the 'near side') whereas the mea- 
surement of jet properties is done for the other parton 
(the 'away side'). Typical examples for measurements 
in which shower biases occur are the near side associate 
hadron distribution for hadron triggered events, or the 
distribution of hadrons inside reconstructed jets. 

In order to study the effect of imposing a trigger con- 
dition on the shower evolution in isolation, we keep the 
parameters which are subject to other biases, i.e. par- 
ton type, initial energy and the strength of the medium 
modification fixed for the following section, i.e. the test 
case is always chosen to be a 20 GeV quark, either in 
vacuum or propagating through a medium such that the 
line-integrated medium-induced virtuality is AQ^ = 5 
GeV 2. 



A. Hard track conditions 

Let us now consider showers in which a trigger condi- 
tion forces a single hard hadron to have the momentum 
Ph- Energy- momentum conservation inside the shower 
requires to recover (in the medium case approximately) 
the original shower-initiating parton energy as E — Pi 
(where hadron masses have been neglected), i.e. in the 
limit where Ph/E is sufficiently large, a significant bias 
for the remaining distribution is found. This is illustrated 
in Fig. [3] where the longitudinal momentum distribution 
inside a vacuum shower is plotted for various values of 
the imposed hard track condition. 
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FIG. 3: Conditional distribution of hadrons at energy E in 
a shower originating from a 20 GeV quark, given a trigger 
hadron with the indicated energy in the same shower. 

In essence, a hard track condition leads (by construc- 
tion) to an enhancement of the hadron yield in the mo- 
mentum region above the track requirement, and by mo- 
mentum conservation to a depletion of the yield below. 
This pattern becomes more pronounced if the trigger 
hadron takes a sizeable fraction of the total jet energy. 
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FIG. 4: Medium over vacuum ratio of conditional hadron 
energy distributions in a shower originating from a 20 GeV 
quark, given a trigger hadron with the indicated energy in the 
same shower. 



turn distribution inside a vacuum shower originating from 
a 20 GeV quark is plotted for various jet energy cuts af- 
ter the shower has been clustered with anti- i^T for the 
indicated radius parameter. 
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Fig. m illustrates how the medium modification of the 
shower responds to a hard track trigger condition at the 
example of the ratio of the longitudinal momentum dis- 
tributions in medium and in vacuum. In the unbiased 
case, a depletion of the yield at high Pt ('jet quenching') 
is balanced by a significant yield increase at low Pt- A 
hard track condition tends to remove the depletion at 
high Pt above the required track momentum. This is a 
very natural outcome — while hard tracks are unlikely in 
the unbiased case and made even less likely by the effect 
of the medium, a single hard track is always guaranteed 
once the trigger condition is imposed, and hence it can- 
not be quenched by the medium. In the presence of such 
a trigger condition, the medium effect may reduce the 
rate of triggered events, but it may no longer lead to a 
quenched high Pt shower pattern. This is a very generic 
finding — imposing a trigger bias always tends to reduce 
the medium modifications of the shower pattern because 
the trigger condition generates protected structures in 
the jet. 



B. Jet energy conditions 

Another commonly found bias on the shower is the 
requirement that a jet with at least a certain energy is 
found. The precise nature of the bias depends on the 
algorithm used to cluster hadrons into jets and their pa- 
rameter settings (often an angular radius parameter R), 
as well as on the energy threshold. Qualitatively, it is 
clear that requiring a substantial flow of the total jet en- 
ergy into a small cone radius selects showers in which 
only few branchings take place and as a consequence the 
longitudinal spectrum is harder than average while the 
transverse shape is more coUimated than average. 

This is shown in Fig. [S] where the longitudinal momen- 



FIG. 5; Conditional distribution of hadrons at energy E in 
a shower originating from a 20 GeV quark, given that the 
shower after clustering with radius R results in a jet energy 
above -Ejet. 

As expected, requiring a very collimated jet by asking 
at least 75% of the jet energy in a cone of radius i? = 0.2 
leads to a sizeable hardening of the hadron spectrum in 
the jet, with the bias successively decreasing for larger 
radii or smaller Ejet- However, unlike in the case of a 
hard track requirements, there are no pronounced dis- 
continuities in the distribution created by a jet energy 
condition. 

20 GeV quark, AQ^ = 5 GeV^ 



• min. bias 




FIG. 6: Medium over vacuum ratio of conditional hadron 
energy distributions in a shower originating from a 20 GeV 
quark, given that the shower after clustering with radius R 
results in a jet energy above Ejet-- 

Fig. m illustrates how the medium modification of the 
shower structure is affected by imposing a jet energy con- 
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dition. As in the case of a hard track condition, gener- 
ically a bias on the shower tends to remove the modifi- 
cation, as an increasingly significant part of the shower 
becomes protected against any modification by the trig- 
ger condition. 



C. Jet mass conditions 

A final state jet property which is currently not used 
as a trigger condition in measurements is the mass of a 
jet. This can be related to the virtuality of the shower- 
initiating parton which in turn determines the phase 
space available for vacuum branchings. Thus, highly vir- 
tual partons undergo a much richer branching history be- 
fore the medium is encountered than hard partons with 
low initial virtuality. In this way, tagging high jet masses 
selects events in which configurations of multiple partons 
undergo medium modification, whereas tagging low mass 
jets tends to prefer configrations which are dominated by 
a single leading parton. The strength of medium modi- 
fications observed in a shower is thus expected to scale 
with the jet mass, which can be exploited to get a more 
differential picture of medium-modified showers [49| . 



VI. CASE STUDY: AWAY SIDE Iaa 

Let us in the following consider a more realistic situa- 
tion in which the trigger condition refers to a pure final 
state condition and hence other types of biases (in this 
section with the exception of the shower bias) occur, in 
particular kinematic, parton type and geometry bias. 

The test case discussed in this section is a measurement 
of Iaa of the conditional yield of hadrons as a function 
of Pt on the away side (which removes the shower bias), 
binned as a function of zt where zt = Pr/Eots- 

Conditional away side yields were first obtained by the 
STAR collaboration^, [13, H^l and the strong quenching 
of the away side correlation peak was almost immediately 
seen as a spectacular confirmation of the expectations of 
monojet events in a medium. Similar measurements have 
now also been performed by the ALICE collaboration at 
LHC kinematics (55| which found somewhat reduced sup- 
pression as compared to the RHIC case. While theoreti- 
cally challenging to compute, dihadron correlations have 
been a valuable tool to probe for instance pathlength- 
dependence of energy loss [H, 113] and to track the fate 
of subleading hadrons [ssj . 

In the case study, -Bobs is always the trigger energy 
given the trigger condition. We consider four different 
trigger trigger conditions: 1) a 7 (7-h) 2) a single hadron 
(h-h) 3) a jet as defined by STAR [l^, including only 
TT"*", 7r~, tt'', K'^, K~,p,p or 7 above 2 GeV clustered with 
a radius of i? = 0.4 (jet-h) and 4) an ideal jet with all 
particles clustered into R — 0.4 (ijet-h). The trigger 
energy range is in all cases 12—15 GeV. 



This selection contains strong kinematical bias (h-h, 
jet-h) as well as weak kinematical bias (ijet-h, 7-h), 
strong parton type bias (7-h, h-h) as well as weak parton 
type bias (jet-h, ijet-h) and strong geometry bias (h-h, 
jet-h) as well as weak geometry bias (7-h, ijet-h). 

There is some freedom in the choice of the away side 
observable, and in principle one could have chosen for 
instance Pr rather than zt- Each of these choices em- 
phasizes different physics: At low Pt, the jet structure 
is determined largely by the appearance of the medium- 
induced radiation. As argued in f50,\ . the enhancement 
region is essentially set by medium physics and thus is 
seen at constant Pt, not constant zt, in which case 
plotting the correlation in Pt emphasizes the relevant 
physics. In contrast, in the high Pt region where z > 0.5, 
energy-momentum conservation is a major influence, and 
the constraints by energy-momentum conservation scale 
on average approximately with zt (i.e. as a constant frac- 
tion of the trigger energy) rather than Pt (i.e. indepdent 
of the trigger energy). The choice made here thus em- 
phasizes the high Pt physics at the expense of obscuring 
the physics of the enhancement due to medium-induced 
radiation. 

Note again that the following results are case studies 
for the sake of illustration rather than model predictions, 
since they do not correct for effects like background fluc- 
tuations in jet finding and ignore the systematic uncer- 
tainty inherent in the choice of the hydrodynamical back- 
ground, which is known to be important in comparison 
with real data 15 li] . 



A. The situation at RHIC 

We consider first the situation for 0-10% central Au-Au 
collisions at RHIC. The distribution of trigger vertices as 
obtained in the model calculation illustrating the amount 
of geometrical bias is shown in Fig. [71 the distribution of 
away side parton momenta given a trigger in the 12-15 
GeV energy range is shown in Fig. [H 

These results confirm in a quantitative way what has 
been stated earlier: Both h-h and jet-h correlations have 
a relatively strong geometry bias to trigger on events in 
which the vertex is close to the surface. This is not so 
for ijet-h correlations (and since the 7 does not undergo 
any final state interaction, the 7-h correlation has no 
geometrical bias at all). 

At the same time, a 7-trigger is, up to intrinsic kT 
smearing, a relatively faithful representation of the par- 
ton kinematics. An ideal jet maps to a somewhat larger 
region in parton pT, whereas jet-h and h-h probe the 
widest range in parton kinematics. 

The parton type bias is summarized in Tab. H] where 
the fraction of gluon jets fgiue is shown on near and away 
side in both vacuum and medium. While the away side 
for the 7-h trigger is almost a pure quark jet sample, 
all other trigger conditions lead to a sizeable gluon jet 
fraction of ^ 60%. 
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FIG. 7: Conditional distribution of production vertices in the transverse plane, given a trigger with observed energy Eobs 
between 12 and 15 GeV in 0-10% central 200 AGeV Au-Au collisions for hadron triggers (left), a jet definition used by STAR 
(middle) and an idealized jet definition (right). In all cases, the trigger object momentum vector defines the —x direction. 



0-10% central 200 AGeV AuAu 

E _^_= 12-15 GeV 
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FIG. 8: Conditional momentum distribution of the away side 
parton given a triggered object in the range of Eobs between 
12 and 15 GeV for various possibilities for the trigger. Shown 
for reference is the situation for p-p collisions (lines) as well 
as the situation in 0-10% central 200 AGeV AuAu collisions 
(symbols). 



Let us briefly review how these biases affect Iaa- A 
strong kinematical bias increases Iaa since the available 
parton energy on the away side increases, giving a larger 
phase space for particle production. Parton type bias 
towards gluon jets decreases Iaa since gluon jets show 
softer fragmentation, a comparison of the numbers sug- 
gests however that in the particular kinematic window 
studied here the differences between vacuum and medium 
are rather small and gluon filtering is not an issue. Fi- 
nally, a strong geometrical bias decreases Iaa since the 
average in-medium pathlength (and hence the strength 
of the medium modifications) grows. However, there is 
no easy a priori argument which would indicate which 
bias determines the end result. 



trigger 


f7ne near 


fZl away / 


near 


/Zf away 


7-h 


N/A 


0.03 


N/A 


0.03 


h-h 


0.04 


0.69 


0.04 


0.69 


jet-h 


0.12 


0.68 


0.08 


0.69 


ijet-h 


0.44 


0.55 


0.33 


0.61 



TABLE I: Conditional fraction of gluon jets on near and away 
side given a trigger object in the range of Eots between 12 and 
15 GeV both in vacuum and in 0-10% central 200 AGeV Au- 
Au collisions. 
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FIG. 9: Away side hadron yield modification as a function of 
zt = Eh/ Eobs for various trigger objects in 0-10% 200 AGeV 
Au-Au collisions. 



The actual outcome of the model in terms of away side 
Iaa is shown in Fig. O Given the fairly sizeable differ- 
ences in geometry and kinematics probed, the outcome 
looks at first glance depressingly similar. All curves show 
suppression at high zt whereas there is enhancement at 
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low zt (which reflects the generic physics of a MMFF as 
determined by comparison with a large body of data [22j 
— energy lost from hard shower modes is recovered in the 
enhanced production of subleading hadrons. Differences 
in detail, in particular between 7-h and ijet-h or h-h and 
jet-h, appear lost. 

There are two possible scenarios which can generate 
such a result: Either a generic effect makes the outcome 
of the computation insensitive to the details of the bias, 
or there is an accidential cancellation of biases acting in 
different directions. 



« y-h default 

•— • y-h, gluon recoil 




FIG. 10: 

The result shown in Fig. [TU] argues that the latter sce- 
nario is true — if the parton type bias is changed to the 
(unphysical) case that only gluons recoil from a 7 trigger, 
the stronger interaction of the gluon with the medium is 
expected to lead to additional softening of the away side 
yield — which is exactly what is observed. Thus, the 
observation that 7-h and ijet-h results fall almost on top 
of each other is not due to some generic mechanism, but 
results from a non-trivial cancellation of biases. 

This in turn argues that if the relative strength of the 
biases can be changed experimentally, the cancellation 
can no longer be expected to occur. One possibility to 
do so is to consider the LHC kinematic range at a signif- 
icantly higher = 2.76 ATeV. 



B. The situation at LHC 

When going from ^ = 200 AGeV to Vs of 2.76 
ATeV with trigger momentum ranges kept fixed, the 
following trends are expected in the biases: First, the 
hard collisions probe the nuclear initial state at lower 



2S„ 



s, and consequently there is a transition to 



a significantly more gluon-dominated regime, as gluons 
increasingly constitute the largest share of the low x par- 
ton distribution. This has an effect on the parton type 
bias. In addition, the momentum spectrum of produced 
partons gets much harder, which implies a weakening of 
the kinematic bias since the 'penalty' for using a very 



energetic parton to produce a high Pt hadron decrases. 
As a result, the correlation between parton momentum 
and jet or leading hadron momentum generically weak- 
ens. Finally, there is also a more copious production of 
bulk matter, both medium temperature and density are 
increasing with y^, which implies a strengthened geomet- 
rical bias. However, since the available kinematic range 
grows ^ whereas the medium density grows as a weak 
power of -^s (for instance ~ -^3"'^^^ in the EKRT model 
[53|V there is some reason to expect a net weakening of 
the geometrical bias. 

Again, note that the following results are for illustra- 
tion and not predictions, as they use a direct extrapola- 
tion from RHIC to LHC energies |52i] with no attempt 
to tune model parameters to LHC data or to explore the 
systematic uncertainty given by the choice of the hydro- 
dynamical background model. 

An explicit computation of the geometrical bias shown 
in Fig. [TT]confirms this expectation — despite the higher 
temperature and density of the LHC medium, the result- 
ing bias on geometry is found to be considerably less due 
to the harder parton spectrum. This can also be seen 
from Fig. [12] where the conditional distribution of away 
side parton momenta given a trigger is shown. 

It is evident that the same range in trigger Pt maps to 
a much wider range in possible parton kinematics at the 
LHC than at RHIC. The underlying reason is again the 
reduced penalty for starting with a high parton energy, 
which in turn is due to the harder primary parton spec- 
trum. The changes in parton type bias are summarized 
in Table HH 



trigger 


r.tue near 


r^tue away / 


"led 

qlue near 


/Xf away 


7-h 


N/A 


0.04 


N/A 


0.04 


h-h 


0.33 


0.79 


0.32 


0.78 


jet-h 


0.47 


0.79 


0.38 


0.80 


ijet-h 


0.77 


0.78 


0.69 


0.78 



TABLE II: Conditional fraction of gluon jets on near and 
away side given a trigger object in the range of Eobs between 
12 and 15 GeV both in vacuum and in 0-10% central 2.76 
ATeV Au-Au collisions. 

As expected, it can be seen that in particularly the 
near side gluon fraction at LHC is much increased over 
RHIC values, but also that the away side gluon fraction 
is more independent of the near side gluon fraction. This 
dependence at RHIC happened because of the dominance 
of the gq — > qg reaction, which is no longer dominating 
at LHC kinematics. 

The final model result in terms of away side Iaa are 
shown in Fig. [131 These results show that the change 
in -^s from RHIC to LHC, resulting in different kine- 
matical, geometrical and parton type biases is in prin- 
ciple strong enough to leave significant traces in observ- 
able quantities. Any similarity between RHIC and LHC 
results should therefore not be seen as caused by the 
same generic (and hence trivial) dynamics, but rather as 
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FIG. 11: Conditional distribution of production vertices in the transverse plane, given a trigger with observed energy Eots 
between 12 and 15 GeV in 0-10% central 2.76 ATeV Pb-Pb collisions for hadron triggers (left), a jet definition used by STAR 
(middle) and an idealized jet definition (right). In all cases, the trigger object momentum vector defines the —x direction. 




FIG. 12: Conditional momentum distribution of the away side 
parton given a triggered object in the range of Eots between 
12 and 15 GeV for various possibilities for the trigger. Shown 
for reference is the situation for p-p collisions (lines) as well 
as the situation in 0-10% central 2.76 ATeV Pb-Pb collisions 
(symbols). Note the change in the scale of the x-axis in com- 
parison with Fig. [S] 



carrying meaningful information in terms of the relative 
strength of biases and their cancellations. 



VII. CASE STUDY — THE LONGITUDINAL 
MOMENTUM DISTRIBUTION OF JETS 

Clustering hadrons into jets has been introduced in 
the study of hard QCD processes in p-p collisions with 
the aim of providing an easy comparison between pQCD 
calculations on the partonic level and the experimentally 
observed hadronic final state. The basic idea is that clus- 
tering largely removes the effect of any soft physics like 



FIG. 13: Away side hadron yield modification as a function of 
ZT — Eh/Eobs for various trigger objects in 0-10% 2.76 ATeV 
Pb-Pb collisions. 



hadronization or additional soft gluon emission which can 
not alter the flux of energy and momentum in the shower 
significantly, and thus a fairly direct comparison of ex- 
perimental observables is made possible. 

It is doubtful if this still constitutes an advantage in the 
study of medium- modified showers, since medium modifi- 
cation occurs predominantly driven by the medium tem- 
perature scale T ^ 300 — 500 MeV (which is soft), i.e. 
clustering into jets tends to suppress the very effect one 
sets out to study. It can be shown that this renders di- 
jet imbalance observations fairly insensitive to even gross 
features of the parton-niedium interaction 58 1. 

One way to overcome this problem is to analyze the 
spectrum of particles in the observed jets and hence get 
a more differential picture. In the language developed 
above, this corresponds to a situation where in addition 
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to kinematic, parton type and geometry bias also the 
shower bias is relevant. The test case considered here is 
an analysis of the longitudinal momentum distribution of 
jets in 2.76 ATeV PbPb collisions clustered from hadrons 
above 1 GeV with anti-^T using R = 0.3 with the jet en- 
ergy Ejet required to fall into the range of 100 - 110 GeV 
(note that this is similar to the fragmentation function 
analysis by the CMS collaboration [SMi)- 

From Fig. [11] we may infer that the geometry bias in 
this situation is weak, and from Fig. [12] we can see that 
we may expect partons from the trigger energy threshold 
Ejet to about 1.5-2 times the trigger energy to contribute 
to the yield of jets in the trigger energy range. The com- 
plication due to the shower bias can be estimated from 
Fig. [5] For parton energies close to Ejet (i-e. parton ener- 
gies around 110 GeV) there is reason to expect a strong 
bias towards a shower structure which is not medium- 
modified, for parton energies sufficiently above Ejet this 
bias gradually lessens, with the relative weight of these 
situations being determined by the combination of kine- 
matical and parton type bias. 
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FIG. 15: Near side hadron yield medium modification in a 
R = 0.3 anti-fcr jet as a function of Pt, shown as full result 
and obtained by neglecting the shower bias. 
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FIG. 14: Conditional momentum distribution of the near side 
parton given a triggered jet in the range of Eobs between 100 
and 110 GeV. Shown for reference is the situation for p-p 
collisions (line) as well as the situation in 0-10% central 2.76 
ATeV Pb-Pb collisions (symbols). 

The kinematic bias as obtained in the model calcula- 
tion is illustrated in Fig. [14] As expected, parton ener- 
gies are probed in a range from about 100 to 150 GeV, 
with a slight shift towards higher energies in the medium 
case. At the same time, the fraction of gluon jets con- 
tributing to the yield in the trigger range decreases from 
fgtue = 0-44 to /g7„^f = 0.36. A large fraction of jets is 
hence required to carry 2/3 of the parton energy inside a 
cone of -R = 0.3, which according to Fig. [5] argues for an 
appreciable shower bias. 

The final result of the model calculation is shown in 
Fig. [15] and compared with a computation in which the 
shower bias effect has been deliberately removed (i.e. the 
fragmentation is computed for a population of shower 



as given by the kinematic, geometry and parton type 
bias as given by evaluating the trigger condition using 
the full simulation, but it is not checked in this run if a 
given shower actually clusters to an Ejet in the trigger 
energy range — note also that without such rejection, 
the available statistics is much higher). 

The results show dramatic differences between taking 
the shower bias into account or not. In all cases, there is 
an enhancement of the yield below Pt ^ 3 GeV (which 
is not properly resolved by the binning) . This is followed 
by a statistically significant region of depletion in the full 
calculation which ends at around 30-40 GeV where the 
full result becomes compatible with unity before statistics 
runs out. Qualitatively, this agrees with CMS measure- 
ments [soj . In contrast, the result without shower bias 
continues to show increasing depletion of the yield up to 
the highest hadron Pt- 

The reason for the peculiar pattern of enhancement, 
depletion and unity observed in the full calculation is a 
good illustration of the interplay between different bi- 
ases. At small Pt, the contribution of gluon jets is still 
appreciable, and so the full calculation shows the same 
enhancement and depletion as the unbiased calculation, 
albeit driven towards unity by the shower bias. How- 
ever, at high Pt the yield is almost exclusively due to 
quark jets since the fragmentation of gluon jets is gener- 
ically softer. Thus, at some point the enhanced fraction 
of quark jets in the medium due to gluon filtering leads 
to a parton type bias towards Iaa > 1 which happens to 
approximately compensate the softening of the spectrum 
due to the medium modification in this kinematic range. 
The net result is Iaa ~ 1 in the high Pt region. 
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VIII. COMPLICATED BIASES 

In several experimentally relevant situations, even 
more complicated bias structures appear. One exam- 
ple are 2+1 triggered correlation in which the trigger 
condition corresponds to observing a coincidence of hard 
hadrons on both the near and the away side. A different 
example are triggered or seeded jets in which clustering 
of the event into jets is only done if a high Pt track 
has been seen in the event. Let us study these cases in 
somewhat more detail. 



A. 2+1 triggered correlations 

While in hadron (or jet) triggered correlations the away 
side parton propagation is constrained in azimuth to be 
approximately back to back with the trigger parton, the 
rapidity of the away side parton is only weakly con- 
strained given the observed rapidity of the near side par- 
ton, and only at sufficiently high Pt kinematics forces 
them to a similar rapidity (see e.g. [g^l)- The original 
motivation for introducing 2+1 triggered correlations in 
which a hard hadron on both the near (Tl) and the away 
side (T2) serves as trigger condition was to explicitly con- 
strain the rapidity of the away side parton, and hence to 
have a better lever arm to study correlations caused by 
energy-deposition into the medium. It was however re- 
alized fairly quickly that such a trigger condition biases 
the event towards minimal medium modification of the 
shower, which tends to make the observation of energy 
redistribution difficult to impossible [Slj . 

Since hard hadron production is rare to begin with, 
the hard fragmentation in coincidence is an even rarer 
phenomenon, and this implies a strong bias in the event 
structure. We may hence expect a strong kinematical 
bias with on average significantly higher parton energies 
probed than in the hadron triggered case, a parton type 
bias leaning towards quark jet coincidences and a sym- 
metric (tangential) geometry bias which minimized the 
in-medium pathlength for both near and away side par- 
ton, combined with a shower bias on each side given by 
the trigger requirement. 

While the original motivation for measuring 2+1 coin- 
cidences has a doubtful prospect of being used in prac- 
tice, 2+1 triggered correlations have the appealing fea- 
ture that changing the momentum range for T2 allows to 
change the underlying bias structure in a profound way 
with minimal effort. The downside is that since hard di- 
hadron coincidences are rare, finite statistics limits their 
usefulness. 

In the following case study, we set Tl to the window 
of 12-15 GeV in order to compare with previous results 
for hadron-triggered correlations and study two ranges 
for T2, 4-8 GeV and 8-10 GeV. We consider the case of 
0-10% central Au-Au collisions at 200 AGeV only. 

Fig. [Tni shows the geometrical bias obtained from the 
model calculation. While in the hadron-triggered case 



there is a surface bias, coming from the requirement of 
having a short in-medium path for the trigger hadron, 
with increased momentum required for T2 this gradually 
changes into a tangential bias for which both near and 
away side parton in-medium pathlength are minimized. 
This means that paths through the center of the medium 
are progressively suppressed and for close to equal mo- 
menta of Tl and T2, chiefly the periphery of the medium 
is probed. 

As expected, the T2 condition has also implications for 
the parton kinematics. This is demonstrated in Fig. 1171 
For the highest T2 range, the mean parton momentum 
probed by the triggered correlation is moved about 10 
GeV higher than for the single hadron trigger. The im- 
plication of this is naturally a substantial suppression of 
the trigger rate. 



trigger /Xe near /„7„", away /^f near /^f away 

h-h 0.04 0.69 0.04 0.69 

T2 = 4-8 GeV 0.071 0.49 0.07 0.38 

T2 = 8-10 GeV 0.10 0.29 0.05 0.20 



TABLE III: Conditional fraction of gluon jets on near and 
away side given a trigger object in the range of Tl between 
12 and 15 GeV and T2 in the indicated range both in vacuum 
and in 0-10% central 200 AGeV Au-Au collisions. 

The evolution of the gluon jet fraction with T2 is shown 
in Tab. IIIII It is apparent that in vacuum the dominance 
of quark jets leading to a near side with correlated gluons 
on the away side is progressively broken. This is a natural 
consequence of the kinematic shift. In the medium, there 
is a strong gluon filtering effect apparent on the away 
side, leading to the dominance of correlated quark jets 
on both near and away side. 

The resulting away side Iaa{zt) for the different 
ranges of T2 is shown in Fig. [THl Perhaps not sur- 
prisingly, a significant enhancement of the yield above 
vacuum is found. This is a result of the strong kine- 
matic bias, shifting parton energy upward and allowing 
for more phase space for subleading hadron production, 
the tangential bias which reduces the medium modifica- 
tion for the away side shower as compared with the sin- 
gle hadron triggered case and the parton type bias which 
drives the away side towards harder quark jets. The im- 
mediate consequence of these biases is a strong reduction 
in the rate at which triggers are produced (which here 
refiects in the larger statistical errors for the dihadron 
triggered results, as significantly more events need to 
be created for this observable than for hadron triggered 
events) . 

B. Jet finding in triggered events 

From an experimental point of view, it is often unde- 
sirable to run jet finding algorithms on a set of minimum 
bias events in heavy-ion collisions, as the vast majority 
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FIG. 16: Conditional distribution of production vertices in the transverse plane, given a dihadron trigger with observed energy 
Tl between 12 and 15 GeV in 0-10% central 200 AGeV Au-Au collisions and T2 set to the indicated values. Shown (left) is 
also the situation without T2 requirement (see Fig. [7|. In all cases, the Tl momentum vector defines the —x direction. 
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FIG. 17: Conditional momentum distribution of the near side 
parton given a dihadron trigger with Tl between 12 and 15 
GeV and T2 in the indicated range. Shown for reference is 
the situation for p-p collisions (line) as well as the situation 
in 0-10% central 200 AGeV Au-Au collisions (symbols). 



of these events will not contain a hard process and hence 
significant numerical effort is used to cluster events which 
are not relevant for the study of hard probes. In such a 
situation, a triggered event sample where events are only 
processed further if they contain a hard track or tower 
(which can be determined early on) can be used. The 
STAR jet analysis [HI exemplifies this strategy for in- 
stance. 

However, triggering on events in this way introduces 
a shower bias. Assumimg that in addition it is required 
that the hard track/tower is part of the leading jet, there 
is a combined bias from both a jet energy and a track 
energy condition. Since the kinematical or geometry bias 
are rather different for jets than for leading hadrons, an 
interesting question is then whether the objects triggered 



0.2 0.4 0.6 0.8 1 

FIG. 18: Away side associate hadron yield medium modifi- 
cation given a single hadron trigger and a dihadron trigger 
with Tl = 12-15 GeV and T2 in the indicated range in 0-10% 
central 200 AGeV Au-Au collisions. 

in this way behave more like jets or like hadrons. 

Obviously this depends on the ratio of clustered jet 
energy to required hadron energy — if the hadron is in 
such a momentum regime that a typical jet contains one 
or more hadrons at this scale, the bias can be expected 
to be small (for instance, requiring a 5 GeV hadron in 
a 100 GeV jet is not expected to bias the jet sample in 
a significant way as the vast majority of 100 GeV jets 
produces hadrons in the 5 GeV range). On the other 
hand, once the hadron carries a significant fraction of the 
total jet energy, jets containing hadrons at such a scale 
become rare and the additional bias will be substantial. 

In Fig. [TOl sample calculations with a combined shower 
bias are shown. The rule of thumb emerging from this 
and similar studies appears to be that the additional bias 
becomes relevant once the hadron energy reaches about 
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20 GeV quark 20 GeV quark 

R = 0.2, E.^i > 15 GeV R = 0.4, E.^^ > 15 GeV 




E [GeV] E [GeV] 

FIG. 19: Conditional distribution of hadrons at energy i5 in a vacuum shower originating from a 20 GeV quark, given that the 
shower clustered using anti-fcr to an energy Ej^t > 15 GeV and a trigger hadron with the indicated energy in the same shower 
for R = 0.2 (left) and R = 0.4 (right). 



half of the jet energy, with a fairly weak dependence on 
the radius used to cluster the jet. Here it has been tacitly 
assumed that kinematics is such that jet finding typically 
recovers 75% of the parton energy, which according to 
Fig. [5] is not a bad assumption. 



IX. CONCLUSIONS 
A. Biases are everywhere 

As the results of this work show, biases occur for al- 
most any high Pt observable that is in any sense related 
to a conditional probability — be it that an explicit trig- 
ger condition is evaluated for the event or be it an implicit 
condition that a jet needs to be clustered before it can be 
analyzed. This means that understanding and discussing 
biases is an integral part of any theoretical analysys of 
hard probes. 

The main structure of the biases involved is usually 
already apparent in the vacuum, and the medium modi- 
fication to the bias structure of the problem can in many 
cases be regarded as a correction. The strength of the 
medium-induced bias is always apparent from the mod- 
ification (in most cases suppression) of the trigger rate, 
which in turn is directly measured in disappearance ob- 
servables such as nuclear modification factors Raa for 
various trigger objects. 

However, the strength of the medium-induced bias 
does not provide an a priori indication of the modifica- 
tion of conditional yield observables — some biases (for 
instance the kinematic bias) may lead to increasing con- 
ditional yields despite a suppression of the trigger rate, 
whereas other biases such as the geometry bias work to- 
wards a suppression of conditional yields. 



B. Biases are important 

As for instance Fig. [T5] indicates, taking biases into ac- 
count properly can change the result of a computation 
quantitatively and even qualitatively. Thus, a theoreti- 
cal model calculation in which the fate of unbiased parton 
showers in a medium is obtained can not be expected to 
compare with data based on the notion that the bias of 
finding a jet is somehow small or would not influence the 
results in a relevant way. In the particular case of the 
longitudinal intra-jet hadron distribution which experi- 
mentally appears unchanged in a medium over a large 
momentum range [s^ , the naive comparison with theory 
without shower bias would find a large discrepancy to the 
data, and hence lead to the conclusion that a previously 
not considered mechanism which makes shower evolution 
in medium similar to vacuum needs to be introduced. 
However, taking the shower bias into account properly, 
the need for any additional mechanism goes away. 

Biases are at least equally important for triggered cor- 
relation measurements — however for these this is usu- 
ally expected, although the relative strength of different 
biases can lead to counter-intuitive results when one e.g. 
expects suppression of a yield based on the geometry bias 
whereas in the actual situation the kinematical bias dom- 
inates, leading to a net enhancement. 



C. Use of biases 

Biases can appear as a nuisance in cases where they 
suppress the physics one is interested in studying, per- 
haps the most striking illustration is the shower bias for 
a hadron or jet trigger (see Figs. S] and ^ where strong 
medium modifications which are a priori present in the 
shower are suppressed by the bias with the effect that Iaa 
is driven towards unity. Such nuisance biases should be 
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avioded if possible — in the context of shower biases, this 
can be done at the simple expense of separating trigger 
side and observable side, i.e. study away side jets with 
a trigger hadron on the near side (as suggested e.g. in 

[H). 

However, in many cases biases can be utihzed by the 
design of a measurement to control the relevant param- 
eters of the hard process to specifically probe the de- 
pendence of a medium-modification on a single control 
parameter. As an example, consider for instance a com- 
parison of jet-h and 7-h correlations at RHIC kinematics. 
According to Fig. [71 the geometrical bias of a sufficiently 
inclusive jet definition is very weak, i.e. in this case 7-h 
and ijet-h correlations probe almost the same geometry. 
According to Fig. [SI they also have a fairly similar kine- 
matical bias, and a small shift in trigger energy range can 
make the underlying parton enery on average the same. 
The main difference between the two situations is then 
given by Tab. [I] from which one can read off that the 7- 
h correlation produces a high fraction of quark jets on 
the away side whereas the ijet-h correlation is dominated 
by away side gluon jets. Thus, in a measurement the 
bias can be designed to specifically probe the different 
evolution of quark and gluon jets in the medium. 

In a similar way, the geometry bias can be systemati- 
cally varied by changing the constituent cut used for the 
clustering (cf. Fig [7]). If the variation of parton kinemat- 
ics associated with the change is compensated for by a 
change in the trigger energy range, the measurement can 
be made to probe various regions of the medium selec- 
tively. 

Of course, such designed observables are inevitably to 
some degree model-dependent. However, in most cases 
the dominant structure of e.g. kinematical or parton 
type bias is given by vacuum QCD, on top of which the 
medium-induced bias is a correction. This means that 



parameters like the necessary shift in trigger energy range 
can be determined approximately by well-known vacuum 
physics and do not have to rely on a particular model of 
parton-medium interaction in a significant way. 



X. SUMMARY 

The ancient Chinese strategist Sun Tzu writes in his 
'Art of War': 

It is said that if you know your enemies and know your- 
self, you will not be imperiled in a hundred battles; if you 
do not know your enemies but do know yourself, you will 
win one and lose one; if you do not know your enemies 
nor yourself, you will be imperiled in every single battle. 

Similarly, one might summarize the results of this work 
as: 

// you understand parton medium interaction and the 
involved biases, everything will become clear. If you un- 
derstand parton-medium interaction but not the biases, 
some observables will make sense, others will appear as 
puzzles; but if you have neither an understanding of bi- 
ases nor a good model of parton-medium interaction, you 
can not know the implication of any hard probe. 
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